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CN ' Abstract 

£ . 

Efficiently computing fast paths in large-scale dynamic road networks (where dynamic traffic 
information is known over a part of the network) is a practical problem faced by several traffic 
information service providers who wish to offer a realistic fast path computation to GPS terminal 
£vq . enabled vehicles. The heuristic solution method we propose is based on a highway hierarchy-based 

shortest path algorithm for static large-scale networks; we maintain a static highway hierarchy and 
perform each query on the dynamically evaluated network. 



O ! 1 Introduction 



Several cars are now fitted with a Global Positioning System (GPS) terminal which gives the exact 
geographic location of the vehicle on the surface of the earth. All of these GPS terminals are now endowed 
with detailed road network databases which allow them to compute the shortest path (in terms of distance) 
between the current vehicle location (source) and another location given by the driver (destination). 
Naturally, drivers are more interested in the source-destination fastest path (i.e. shortest in terms of 
■ travelling time). The greatest difficulty to overcome is that the travelling time depends heavily on the 

amount of traffic on the chosen road. Currently, some state agencies as well as commercial enterprises 
are charged with monitoring the traffic situation in certain pre-determined strategic places. Furthermore, 
traffic reports are collected from police cars as well as some taxi services. The dynamic traffic information, 
however, is as yet limited to a small proportion of the whole road network. 



The problem faced by traffic information providers is currently that of offering GPS terminal enabled 
drivers a source-destination path subject to the following constraints: (a) the path should be fast in terms 
of travelling time subject to dynamic traffic information being available on part of the road network; (b) 
traffic information data are updated approximately each minute; (c) answers to path queries should be 
computed in real time. Given the data communication time and other overheads, constraint (c) practically 
asks for a shortest path computation time of no more than 1 second. Constraint (b) poses a serious 
problem, because it implies that the fastest source-destination path may change each minute, giving 
an on-line dimension to the problem. A source-destination query spanning several hundred kilometers, 
which would take several hours to travel, would need a system recomputing the fastest path each minute; 
this in turn would mean keeping track of each query for potentially several hours. As the estimated 
computational cost of this requirement is superior to the resources usually devoted to the task, a system 
based on dynamic traffic information will not, in practice, ever compute the on-line fastest path. As a 
typical national road network for a large European country usually counts several million junctions and 
road segments, constraint (c) implies that a straight Dijkstra's algorithm is not a viable option. In view of 
constraint (a), in our solution method fast paths can be efficiently computed by means of a point-to-point 
hierarchy-based shortest path algorithm for static large-scale networks, where the hierarchy is built using 
static information and each query is answered on the dynamically evaluated network. 
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This paper makes two original scientific contributions (i) We extend a known hierarchy-based shortest 
path algorithm for static large-scale undirected graphs (the Highway Hierarchies algorithm |SS05j ) to 
the directed case. The method has been developed and tested on real road network data taken from 
the TeleAtlas France database |NV05j . We note that the original authors of [SS05j have extended the 
algorithm to work on directed graphs in a slightly different way than ours (see [SS06]). (ii) We propose a 
method for efficiently finding fast paths on a large-scale dynamic road network where arc travelling times 
are updated in quasi real-time (meaning very often but not continuously). 

In the rest of this section, we discuss the state of the art as regards shortest path algorithms in dynamic 
and large-scale networks, and we describe the proposed solution. The rest of the paper is organized as 
follows. In Scction[2]we briefly review the highway hierarchy-based shortest path algorithm for static large- 
scale networks, which is one of the important building blocks of our method, and discuss the extension of 
the existing shortest-path algorithm to the directed case. Section [3] discusses the computational results, 
and Section U] concludes the paper. 

1.1 Shortest path algorithms in road networks 

The problem of computing fastest paths in graphs whose arc weights change over time is termed the 
Dynamic Shortest Path Problem (DSPP) [BRTedj . The work that laid the foundations for solving 
the DSPP is |CH66| (a good review of this paper can be found in |Dre69j . p. 407): Dijkstra's algorithm 
is extended to the dynamic case through a recursion formula based on the assumption that the network 
G = (V, A) has the FIFO property: for each pair of time instants t, t' with t < t': 

V (u,v) e A r uv (t) + t < T uv (t') + 1', 

where T uv (t) is the travelling time on the arc (u,v) starting from u at time t. The FIFO property 
is also called the non- overtaking property, because it basically says that if A leaves u at time t and 
B at time t' > t, B cannot arrive at v before A using the arc (u,v). The shortest path problem in 
dynamic FIFO networks is therefore polynomially solvable [Cha98 , even in the presence of traffic lights 
AOPS03 . Dijkstra's algorithm applied to dynamic FIFO networks has been optimized in various ways 
BRTcd, Cha98]; the A* one-to-one shortest path algorithm has also been extended to dynamic networks 
[CS02] , The DSPP is NP-hard in non-FIFO networks |Dea04j . 

Although in this paper we do not assume any knowledge about the statistical distribution of the arc 
weights in time, it is worth mentioning that a considerable amount of work has been carried out for 
computing shortest paths in stochastic networks. A good review is jFHK+05] , 

The computation of exact shortest paths in large-scale static networks has received a good deal of 
attention CZOl]. The established practice is to delegate a considerable amount of computation to a 
preprocessing phase (which may be very slow) and then perform fast source-destination shortest path 
queries on the pre-processed data. Recently, the concept of highway hierarchy was proposed in [SS051 
ISch05[ ISS06] . A highway hierarchy of L levels of a graph G = (V, A) is a sequence of graphs G = 
G°,...,G L with vertex sets V° = V, V 1 D . . . D V L and arc sets A = A, A 1 D . . . D A L ; each arc has 
maximum hierarchy level (the maximum i such that it belongs to A 1 ) such that for all pairs of vertices 
there exists between them a shortest path (ai, . . . , <Xfc), where a,i are the consecutive path arcs, whose 
search level first increases and then decreases, and each arc's search level is not greater than its maximum 
hierarchy level. A more precise description is given in Section^ The A* algorithm has also been extended 
to use a concept, reach, which has turned out to be closely related to highway hierarchies (see IGKW05] ). 

1.2 Description of the solution method 

The solution method we propose in this paper efficiently finds fast paths by deploying Dijkstra-like 
queries on a highway hierarchy built using the static arc weights found in the road network database, 
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but used with the dynamic arc weights reflecting quasi real-time traffic observations. This implies using 
two main building blocks: highway hierarchy construction (the Highway Hierarchie^] algorithm extended 
to directed graphs), and the query algorithm. Consequently, the implementation is a complex piece of 
software whose architecture has been detailed in the appendix. 

• Highway hierarchy. Apply the directed graph extension of the HH algorithm (sec Section [2|) to 
construct a highway hierarchy using the static road network information. In particular, arc travelling 
times are average estimations found in the database. This is a preprocessing step that has to be 
performed only when the topology of the road network changes. The CPU time taken for this step 
is not an issue. 

• Efficient path queries. Efficiently address source-destination fast path requests by employing a 
multi-level bidirectional Dijkstra's algorithm on the dynamically evaluated graph using the highway 
hierarchy structure constructed during preprocessing. This algorithm is carried out each time a path 
request is issued; its running time must be as fast as possible, in any case not over 1 second. 



2 Highway Hierarchies algorithm on dynamic directed graphs 

The Highway Hierarchies algorithm [SS05, Sch05] is a fast, hierarchy-based shortest paths algorithm 
which works on static undirected graphs. HH algorithm is specially suited to efficiently finding shortest 
paths in large-scale networks. Since the HH algorithm is one of our main building blocks, we briefly 
review the necessary concepts. 

The Highway Hierarchies algorithm is heavily based on Dijkstra's algorithm |Dij59| , which finds the 
tree of all shortest paths from a root vertex r to all other vertices v s V of a weighted digraph G = (V, A) 
by maintaining a heap H of reached vertices u with their associated (current) shortest path length c(u) 
(elements of the heap are denoted by [u, c(u)]. Vertices which have not yet entered the heap (i.e. which are 
still unvisited) are unreached, and vertices which have already exited the heap (i.e. for which a shortest 
path has already been found) are settled. Dijkstra's algorithm is as follows. 

1. Let H = {[r,0]}. 

2. If H = 0, terminate. 

3. Let u be the vertex in H with minimum associated path length c(u). 

4. Let H = H \ {u}. 

5. For all v E 5 + (u), if c(u) + t uv < c(v) then let H = (H \ {[v, c(«)]}) U {[«, c{u) + t uv ]}. 

6. Go to 2. 

A bidirectional Dijkstra algorithm works by keeping track of two Dijkstra search scopes: one from the 
source, and one from the destination working on the reverse graph. When the two search scopes meet 
it can be shown that the shortest path passes through a vertex that has been reached from both nodes 
( [Sch05j . p. 30). A set of shortest paths is canonical] if, for any shortest pathp = (u\, . . . , Uj, . . . , Uj . . . , Uk) 
in the set, the canonical shortest path between Ui and Uj is a subpath of p. 

The HH algorithm works in two stages: a time-consuming pre-processing stage to be carried out only 
once, and a fast query stage to be executed at each shortest path request. Let G° — G. During the first 
stage, a highway hierarchy is constructed, where each hierarchy level G , for / < L, is a modified subgraph 

^^From now on, simply HH 

2 Dijkstra's algorithm can easily be modified to output a canonical shortest paths tree (see ISch0 5 . Appendix A.l — can 
be downloaded from http://algo2.iti.uka.de/schultes/hwy/ ). 
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of the previous level graph G 1 ^ 1 such that no canonical shortest path in G i_1 lies entirely outside the 
current level for all sufficiently distant path endpoints: this ensures that all queries between far endpoints 
on level / — 1 are mostly carried out on level I, which is smaller, thus speeding up the search. Each shortest 
path query is executed by a multi-level bidirectional Dijkstra algorithm: two searches are started from the 
source and from the destination, and the query is completed shortly after the search scopes have met; at 
no time do the search scopes decrease hierarchical level. Intuitively, path optimality is due to the fact that 
by hierarchy construction there exist no canonical shortest path of the form (ax, . . . , a{, . . . , aj . . . , dfe . . .), 
where aj, Oj, au S A and the search level of aj is lower than the level of both aj, dk] besides, each arc's 
search level is always lower or equal to that arc's maximum level, which is computed during the hierarchy 
construction phase and is equal to the maximum level I such that the arc belongs to G l . The speed of 
the query is due to the fact that the search scopes occur mostly on a high hierarchy level, with fewer arcs 
and nodes than in the original graph. 

2.1 Highway hierarchy 

As the construction of the highway hierarchy is the most complicated part of HH algorithm, we endeavour 
to explain its main traits in more detail. Given a local extensionality parameter H (which measures 
the degree at which shortest path queries are satisfied without stepping up hierarchical levels) and the 
maximum number of hierarchy levels L, the iterative method to build the next highway level / + 1 starting 
from a given level graph G l is as follows: 

1. For each v £ V, build the neighbourhood N l H (v) of all vertices reached from v with a simple Dijkstra 
search in the 2-th level graph up to and including the H-st settled vertex. This defines the local 
extensionality of each vertex, i.e. the extent to which the query "stays on level I" . 

2. For each v G V: 

(a) Build a partial shortest path tree B(v) from v, assigning a status to each vertex. The initial 
status for v is "active" . The vertex status is inherited from the parent vertex whenever a 
vertex is reached or settled. A vertex w which is settled on the shortest path (v, u, . . . , w) 
(where v ^ u ^ w) becomes "passive" if 

\N l H (u)nN l H (w)\<l. (1) 

The partial shortest path tree is complete when there are no more active reached but unsettled 
vertices left. 

(b) From each leaf t of B(v), iterate backwards along the branch from t to v. all arcs (u, w) such 
that u g" N l H {t) and w g N l H (v), as well as their adjacent vertices u,w, are raised to the next 
hierarchy level l + l. 

3. Select a set of bypassable nodes on level l + l; intuitively, these nodes have low degree, so that the 
benefit of skipping them during a search outweights the drawbacks (i.e., the fact that we have to 
add shortcuts to preserve the algorithm's correctness). Specifically, for a given set C Vi+i 
of bypassable nodes, we define the set Si+i of shortcut edges that bypass the nodes in Bj+i: for 
each path p = (s, bi, b%, ■ ■ ■ , bk,t) with s,t € VJ+x \ Bl and 6j € 1 < i < k, the set Si+\ 
contains an edge (s,t) with c(s,t) — c(p). The core G' l+1 = (V l ' +l7 E' l+1 ) of level I + 1 is defined 
as:^' +1 = V l+1 s B l+ i, E' l+1 = (E l+1 n (V/ +1 x Vj' +1 )) U S l+1 . 

The result of the contraction is the contracted highway network G\ , l5 which can be used as input for 
the following iteration of the construction procedure. It is worth noting that higher level graphs may be 
disconnected even though the original graph is connected. 

2.1 Example 

Take the directed graph G — (V, A) given in Fig. [I] (above). We are going to construct a road hierarchy 
with H = 3 and L = 1 on G. First we compute N°(v) for all v £ V — {vq, . . . , v§}. 
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v 



N$(v) 



Vo,Vl,V 2 



{0,1,2} 
{3,4,5} 
{3,5,6} 



Next, we compute B(v) for all v e V and raise the hierarchy level of the relevant arcs from the leaves to 
B(v) to v. We only discuss the computation of B(vq) in detail as the others are similar. 

1. Vertex vq is initialized as an active vertex. 

2. Dijkstra's algorithm is started. 

(a) vq is settled (cost 0) on the empty path, so the passivity condition |ip does not apply; 

(b) v\ and v 2 are reached from vq with costs resp. 1 and 2, and inherit its active status; 

(c) vi is settled (cost 1) on the path (vq, v\) and condition |ip does not apply; 

(d) vq is reached from v\ with cost 1 + 4 = 5 and set to active; 

(e) v 2 (cost 2) is settled on (vq,v 2 ); 

(f) V4 is reached from 2 with cost 2 + 6 = 8 and set to active; 

(g) Vq (cost 5) is settled on the path (v ,vi,vg): since N°(vi) f] N°(vq) — 0, condition |Ip is 
verified, and v§ is labeled passive; 

(h) V3 is reached from v§ with cost 1 + 4 + 4 = 9 and set to passive. 

(i) V4 (cost 8) is settled on the path (v ,v 2 ,V4): since N®(v 2 ) n N®(v4,) = 0, condition |Ip is 
verified, and V4 is labeled passive; 

(j) V5 is reached from V4 with cost 2 + 6 + 2 = 10 and set to passive; 

(k) the only unsettled vertices are v 3 and v 5 . Since both are reached and passive, the search 
terminates. 

3. The leaf vertices of B(vo) are V4 and vq. 

(a) From t — V4, we iterate backwards along the arcs on the path (vq, v 2 , V4): the arc (v 2 , V4) has 
the property that v 2 ^ N§(vi) and V4 $ N®(v 2 ), so its hierarchy level is raised to I + 1 = 1 
(the other arc on the path, (vq,v 2 ), stays at level 1 = 0); 

(b) from t = vq, we iterate backwards along the arcs on the path (vo,vi,Ve): the arc (ui,Ug) -has 
the property that v\ £ N®(v§) and vq ^ N®(vi), so its hierarchy level is raised to 1 (the other 
arc on the path stays at level ). 

Fig.^shows the hierarchy at level 1. 

2.2 Extension to directed graphs 

The original description of the HH algorithm SS05] applies to undirected graphs only; in this section we 
provide an extension to the directed case. It should be noted that the HH algorithm was extended to the 
directed case by the authors (see |SS06| ) in a way which is very similar to that described here. However, 
we believe our slightly different exposition helps to clarify these ideas considerably 

The algorithm for hierarchy construction, as explained in Section I2TT1 works with both undirected and 
directed graphs. However, storing all neighbourhoods N l H (v) for each v and I has prohibitive memory 
requirements. Thus, the original HH implementation for checking whether a vertex v is in N l H (u) is 
based on comparing the distance d(u, v) with the "distance-to-border" (also called slack) from u to the 
border of its neighbourhood N l H (u). The "distance-to-border" d l H (u) is a measure of a neighbourhood's 
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Figure 1: The graph of Example 1 2. II (left) and its highway hierarchy for H = 3, L = 1 (right): the dashed 
lines indicate arcs at level 0, the solid lines indicate arcs at level 1. 

radius, and is defined as the distance d(u,v) where v is the farthest node in N l H (u), i.e. the cost of the 
shortest path from u to the H-th settled node when applying Dijkstra's algorithm on node u at level I. 
This is the basis of the slack-based method in [SchOS . p. 19 (from which we draw our notation). In the 
partial shortest paths tree B(sq) computed in Step[5a|of the algorithm in Section |2~TI the slack A(u) is 
recursively computed for all u £ B(sq) starting from the leaves to of B(sq), as follows. 

1. Initialise a FIFO queue Q to contain all nodes u of B(sq), ordering them from the farthest one to 
the nearest one with respect to sq. 

2. Set A (it) = d l H {u) for u a leaf of B(sq) and +oo otherwise. 

3. If Q is empty, terminate. 

4. Remove u from Q, and let p be its predecessor in B(so)- 

5. If A(p) = +oo and p £" N l H (so), p is added to Q. 

6. Let A(p) = min(A(p), A(u) - d(p, u)). 

7. If A(p) < 0, the edge (p,u) is lifted to the higher hierarchical level. 

8. Return to Step OH 



The algorithm works because Thm. 2 in |Sch05j proves that condition A(p) < is equivalent to the 
condition of Step [2b] of the algorithm in Section 12.11 The cited theorem is based on the following 
assumption: 

Vu£V (u # N l H (t ) -» d l H (t ) - d(u, t ) < 0). (2) 
This condition may fail to hold for directed graphs, since d(u, to) ^ d(to, u). 

To make Assumption [2] hold, we have to consider a neighbourhood radius computed on the reverse 
graph, that is the graph G = (V,A) such that (u,v) £ A (v,u) £ A. Thus, we modified the original 

implementation to compute, for each node, a reverse neighbourhood N H (v) (see Figure [2]), so that we 

are able to store the corresponding reverse neighbourhood radius d H (u)^u £ V. We replace Step-in 
the algorithm above by: 



[2ji. Set A(u) = d H (u) for u a leaf of -B(so) an d +oo otherwise. 



We are now going to prove our key lemma. 
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2.2 Lemma 

Let u,s eV and t a leaf in B(s). If u £ N H (t) then d H (t) - d(u, t) < 0. 

Proof. Suppose d(u,t) < cl H {t). By definition, this means that there is a shortest path in N l H {t) which 
connects u to t. Therefore, u S N l H {t) against the hypothesis. □ 

It is now straightforward to amend Thm. 2 in [Sch05j to hold in the directed case; all other theorems 

in |Sch05| need similar modifications, replacing N l H (t) with N l H (t) and d l H {t) with d H (t) whenever t is 
target node or is "on the right side" of a path - it will always be clear from the context. The query 
algorithm must me modified to cope with these differences, using d H (t) instead of d l H (t) whenever we are 
searching in the backwards direction. 




Figure 2: An example which shows neighbourhoods and reverse neighbourhoods with H — 3; only solid 
arcs are lifted to a higher level in the hierarchy. Note that arcs (p,t) and (p',t) are not lifted even if 

p,p' ^ N l H (t); this is because p,p' £ N l H (t), and for target node we consider the reverse neighbourhood. 

Interestingly, the problem with the slack-based method was first detected when our original imple- 
mentation of the HH algorithm failed to construct a correct hierarchy for the Paris urban area. This 
shows that the extension of the algorithm to the directed case actually arises from real needs. 

2.3 Heuristic application to dynamic networks 

The original Highway Hierarchies algorithm, as described above, finds shortest paths in networks whose 
arc weights do not change in time. By forsaking the optimality guarantee, we adapt the algorithm to the 
case of networks whose arc weights are updated in quasi real-time. Whereas the highway hierarchy is 
constructed using the static arc travelling times from the road network database, each point-to-point path 
query is deployed on a dynamically evaluated version of the highway hierarchy where the arcs are weighted 
using the quasi real-time traffic information. In particular, in all tests that involved a comparison with 
neighbourhood radius we use the static arc travelling times, while for all evaluations of path lengths or 
of node distances we use the real-time (dynamic) travelling times. This means that the static travelling 
times are used to determine neighbourhood's crossings, and thus to determine when to switch to a higher 
level in the hierarchy, while the key for the priority queue for HH algorithm is computed using only 
dynamic travelling times. 



3 Computational results 

In this section we discuss the computational results obtained by our implementation. As there seems 
to be no other readily available software with equivalent functionality, the computational results are not 
comparative. However, we establish the quality of the heuristic solutions by comparing them against 



3 COMPUTATIONAL RESULTS 



8 



the fastest paths found by a plain Dijkstra's algorithm. We mention here two different approaches: 
dynamic highway- node routing ([SS07J, which uses a selection of nodes operated by the HH algorithm to 
build an overlay graph (see HSW06J), and dynamic ALT ([DW07 ), which is a dynamic landmark-based 
implementation of A* . Both approaches, however, although very performing with respect to query times, 
require a computationally heavy update phase (which takes time in the order of minutes), and thus are 
not suitable for our scenario, where, supposedly, each arc can have its cost changed every 2 minutes 
(roughly). 

We performed the tests on the entire road network of France, using a highway hierarchy with H = 65 
and L = 9. The original network has 7778913 junctions and 17154042 road segments; the number of 
nodes and arcs in each level is as follows. 



level 





1 


2 


3 


4 


5 


6 


7 


8 


9 


nodes 


7778913 


1517291 


433286 


182474 


91888 


53376 


34116 


23338 


16445 


11790 


arcs 


17154042 


3461385 


1283000 


583380 


308249 


183659 


119524 


81170 


57235 


41092 



We show the results for queries on the full graph without dynamic travelling times in Table [TJ in 
this case, all paths computed with the HH algorithm are fastest paths. In Table [2j instead, we record 
our results on a graph with dynamic travelling times; we also report the relative distance of the solution 
found with our heuristic version of the HH algorithm and the fastest path computed with Dijkstra, and, 
for comparative reasons, the results of the naive approach which consists in computing the traffic-free 
optimal solution with the HH algorithm (i.e., on the static graph) and then applying dynamic times on 
the so-found solution. Dynamic travelling times were taken choosing, for each query, one out of five sets 
of values recorded in different times of the day for each of the 29384 arcs with dynamic information. 

Although this number is small with respect to the total number of arcs in the graph, it should be 
noted that most of these arcs correspond to very important road segments (highways and national roads). 
All arcs that did not have a dynamic travelling time were assigned a different weight at each query, 
chosen at random with a uniform distribution over [ry, 15ry ], where Ty is the reference time for arc 
This choice has been made in order to recreate a difficult scenario for the query algorithm: even if the 
number of arcs with real traffic information is still small, it is going to increase rapidly as the means for 
obtaining dynamic information increase (e.g. number of road cameras, etc.), and thus, to simulate an 
instance where most arcs have their travelling time changed several times per hour, we generated each 
arc's cost at random. The interval [r^, 15Tjj] is simply a rough estimation of a likely cost interval, based 
on the analysis of historical data. All tables report average values over 5000 queries. All computational 
results in Table [T] and [2] have been obtained on a multiprocessor Intel Xeon 2.6 GHz with 8GB RAM 
running Microsoft Windows Server 2003, compiling with Miscrosoft Visual Studio 2005 and optimization 
level 2. 

Computational results show that, although with no guarantee of optimality, our heuristic version 
of the algorithm works well in practice, with 0.55% average deviation from the optimal solution and a 
recorded maximum deviation of 17.59%; query times do not seem to be influenced by our changes with 
respect to the original version of the algorithm. The naive approach of computing the shortest path 
on the static graph, and then applying dynamic times, records an average error of 2%, but it has a 
much higher variance, and a maximum error of 27.95%; although the average error is not high, it's still 
almost 4 times the average error of the more sofisticated approach, and the high variance suggests lack of 
stability in the solution's quality. The low value recorded for the average error with the naive approach (in 
absolute terms) can be explained as a consequence of the following two facts: travelling times generated 
at random on arcs without real-time traffic information cannot simulate real traffic situation, because 
they lack spatial coherence (i.e. they do not simulate congested nearby zones) and traffic behaviour 
information (i.e. the fact that during peak hours important road segments are likely to be congested, 
while less important roads are not), thus making the task of finding a fast path easier; besides, the 
average query on such a large graph corresponds to a very long path (296 minutes on the traffic-free 
graph, 2356 minutes on the dynamic graph), and on long paths it is usually necessary to use highways 
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or national roads regardless of their congestion status - which is exactly what the HH algorithm does. 
This last sentence is supported by the fact that, if we consider only the 500 shortest queries in terms of 
path length, the average error of the naive approach increases to 3.60%, while the average error of the 
heuristic version increases to 0.97%; this is in accord with the fact that on short paths the influence of 
traffic is greater, because alternative routes that do not use highways are more appealing, while on long 
paths using highways is often a necessary step. However, in relative terms, the heuristic version of the 
HH algorithm performs significantly better than the naive approach proposed for comparison, and we 
expect the difference to increase (in favour of the heuristic algorithm) if applied to a graph fully covered 
with real traffic information. 

Figure [3] shows how the optimal and the heuristic path may differ; since the hierarchy built on the 
static graph emphasizes important roads, the heuristic algorithm applied on the dynamically weighted 
graph still tends to use highways and national roads even when they are congested (up to a certain 
degree), thus sometimes losing optimality 





Dijkstra's algorithm 


HH algorithm 


# settled nodes 


2275563 


18966 


# explored nodes 


2587112 


36200 


query time [sec] 


11.830 


0.099 



Table 1: Computational results on the static graph: average values 





Dijkstra's algorithm 


HH algorithm 

naive approach 


HH algorithm 

heuristic version 


# settled nodes 


2280872 


19174 


19099 


# explored nodes 


2594361 


36581 


36492 


query time [sec] 


11.917 


0.100 


0.099 


distance from optimum (variance) 


0% 


2.00% (5.00) 


0.55% (0.45) 



Table 2: Computational results on the graph with dynamic times: average values 



4 Conclusion 

We present a heuristic algorithm for efficiently finding fast paths in large-scale partially dynamically 
weighted road networks, and benchmark its application on real-world data. The proposed solution is 
based on fast multi-level bidirectional Dijkstra queries on a highway hierarchy built on the statically 
weighted version of the network using the Highway Hierarchies algorithm, and deployed using the dynamic 
arc weights. Computational results show that, although with no guarantee of optimality, the proposed 
solution works well in practice, computing near-optimal fast paths quickly enough for our purposes. 
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Figure 3: Fast paths calculated with different algorithms; each number identifies a path, paths are par- 
tially overlapping. 1: Dijkstra's algorithm (optimal solution) with real-time arc costs; dynamic travelling 
time: 24 minutes and 6 seconds. 2: HH algorithm (heuristic solution) with real-time arc costs; dynamic 
travelling time: 25 minutes and 5 seconds. 3: HH algorithm without real-time arc costs (traffic-free 
optimal solution); dynamic travelling time: 37 minutes and 5 seconds. 
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